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Abstract 

We argue that the large- N c behaviour of the Yukawa couplings in the 
Skyrme model involves issues more subtle than the vanishing of linear 
fluctuations needed for classical stability of the skyrmion. The chiral 
fluctuations about the skyrmion must be quantized in order to reach a 
conclusion. An improved quantization procedure allows us to confront 
this question directly 

The pion-nucleon coupling constants g n NN (CP conserving) and g nNN 
(CP violating) are calculated in the large- N c , three-flavour Skyrme model 
by direct evaluation of the leading matrix elements appearing in the 

3 

LSZ reduction formula. We find that g n NN ~ Nc , but that, at most, 
_ i 

~9irNN ~ iti^Nc 2 . These results show that the leading contribution to the 
neutron electric dipole moment in large- N c Skyrme model is the direct 
one (D n ~ N c m^), rather than the pion loop contribution. 



* Supported in part by the DOE under grant DE-FG02-92ER40706 



1. Introduction 



Some time asp one of us considered chiral fluctuations about a rotating, three- 
flavour skyrmiona 

U(r,t) = U+(T,t) A(t) U (r) A~\t) U^t) (1.1) 

where U (r) is a static solution of the U(3) Skyrme model, A(t) is an SU(3) matrix 
of collective coordinates, and the matrix field (summation convention throughout) 

Uj, = exp{— 4> a \ a } (1.2) 

is built from the nonet of pseudoscalar meson fields <p a (a = 0, 1, . . . , 8) and the 



Gell-Mann matrices A a (a = 1, ... ,8) with A = y 2/31. Insertion of expression [LI 
into the Skyrme lagrangian, £ s , yields the elegant result 

C s (U(t, £)) = CsiAUoA- 1 ) + C a (U 24> ) + Cj^AUoA- 1 - 0) , (1.3) 

in which the right-hand side is isomorphic to the low-energy soft-meson and baryon 
effective lagrangian 

C eS = C B (tlj) + C a (n) + C I (^,n) (1.4) 

which describes free baryons interacting with the mesons of the non-linear cr-model. 
Since |1.4| obeys all the soft-meson theorems of chiral symmetry, one expects that 
these theorems are obeyed by the skyrmion- meson system of |1.3| as well. 

The Skyrme lagrangian gives a formal result for the Yukawa couplings of the 
skyrmion-meson interaction. However, to arrive at a conclusive result for the Yukawa 
couplings involves a number of subtleties which lead to the central issues of this pa- 



per. Given that |L3] includes appropriate CP violating terms, one can compute both 
the CP violating pion-nucleon scalar coupling constant g nNN and the CP conserving 
pseudoscalar coupling constant q^nn appearing in the interaction term of eq. Ol 



d(ip; tt) = 7T • ipr(i"f 5 g nNN + g nNN )ip . (1.5) 

Since the Skyrme model is formulated in terms of both chiral symmetry and the 
large N c (number of colors) expansion, one expects that its predictions for q^nn 
and g nNN are consistent with QCD in the same limits. Naively, the dependence of 
the couplings on the pion (nucleon) mass m n (Mn) and the number of colors N c is, 
from |1.3| and ref . [I], 

g^NN ~ n! g nNN ~ mlN? M N ~N C . (1.6) 

Current algebra provides an ostensibly model-independent derivation of the leading 
contribution to the electric dipole moment of the neutron D n in the chiral limit! 

n g-wNN g^NN i ,M N 

D n = . 2B/f . In . 1.7 
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If eq. [T]6| were correct, the leading contribution in the combined chiral and large N c 
limit (and so dominant prediction of the Skyrme model in these limits) would be 
D n ~ N c ml lnm n . 

Recently significant progress in resolving this issue has been made in calculating 
g-irNN (in the two flavour Skyrme model) by the authors of ref. 0. It is the purpose 
of this paper to calculate the CP conserving and and CP violating 7r-N couplings 
9-kNN, and g nNN by the methods of ref. generalized to the three-flavour Skyrme 
model. We will show that although the pseudoscalar coupling g n NN indeed has the 



leading behaviour of eq. O, the CP violating coupling g nNN is down by a factor of 



at least N c 1 from that of eq. |1.6| . Therefore, the pion contribution to D n is at most 

Dn^milnm^ (1.8) 



and is therefore subdominant to the direct CP violating behaviour D n ~ m^Nc 
computed in ref. [|. (For a review see ref. |5|.) 

It is important to understand why there is an issue with regard to the Skyrme 
model Yukawa couplings. These couplings embody the linear fluctuations about 
the Skyrme classical solution. For constant collective coordinates (A = constant), 
this linear coupling must vanish for the classical solution to be stable under linear 
fluctuations. Thus the classical, linear meson-skyrmion interactions behave as 

£~<p{A terms) (1.9) 

if the Skyrme equations of motion are imposed. The time derivatives in eq. Ol imply 



a suppression of at least a factor of N~ l compared to the behaviour given in |1.6| . If 
the chiral fluctuations are only considered classically, this would be the end of the 
story. However, the classical argument cannot be trusted since the chiral fluctuations 
must be quantized in order to properly account for the mesons as quantum fields.!' i 
This is the origin of a number of complications. 

The canonical momenta of the rotating skyrmion and associated mesons are 
not independent, and the separation into components transverse to the zero-mode 
solutions is "gauge dependent." Correctly calculating an S'-matrix element — which 
must be gauge invariant — such as the 7r-N scattering amplitude is a complicated 
task.i The advance made by the authors of ref. |3] is to calculate the Born terms of 
7T-N scattering in a manifestly gauge-invariant manner in the (two-flavour) Skyrme 
model. They obtain the same result as that given by the gauge-dependent separation 
of skyrmion and meson momenta,i confirming that the Yukawa coupling g^NN has 
the large N c behaviour of eq. [Tl| This suggests that the large- N c analysis of CP 
conserving 7T-N scattering in the Skyrme model is consistent with the implications 
of chiral symmetry. 

The gauge-invariant calculation presented in ref. |3| is considerably simpler than 
the calculations in a gauge that separates the skyrmion and meson momenta. As 
a result, application of this method to the calculation of the tt — N couplings in 
the three-flavour Skyrme model yields a straightforward estimate of the large-iV c 
behaviour of g n NN and ~g TNN . 
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2. The CP Violating Skyrmion 

In order to obtain strong CP violation, one must consider a U(3) non-linear er- 
model so as to include (i) an effective interaction representing the U(l) anomaly 
and (ii) a CP violating mass- matrix. The appropriate effective lagrangian in the 
large N c limit is then! 

C{U) = C S (U) + ^Tt(MU + -M- Mt) + -^(Tr[ln U - In U ] ]f + C wz 

lb o4iv c 

(2.1) 

where Cwz is the Wess-Zumino 5-form, U is a U(3) matrix field, a depends on the 
physical meson masses, 

£s(V) = %^mU){d^)\ + ^e-'Tr{[(d,U)W , {d u U)U^} (2.2) 

and 

M = diag[M 1 exp(z^-),M 2 exp(^^),M 3 exp(ze^)] . (2.3) 

In order to prevent strong CP violation from creating Goldstone bosons from the 
vacuum we require!' @ 

ft = ft = & = 

a + M 3 ( 2 - 4 ) 
6 a + Mi 

up to terms of 0((3 2 ). It is convenient to define 9 = arg detM, the effective CP 
angle, and to set 

e = x °i-k + k + w) (2 - 5) 

so that Xg = (3. 

The classical SU(3) Skyrme ansatz, neglecting CP violation, is00 

tt i \ ( c °s F + it ■ r sin F „ | \/2, „ . 1 -W . . 
t/b(r) = = — (cosF-l)A 8 +-(2cosF+l)/+zr-rsmF 

V oo i J 3 3 

(2-6) 

where F = F(r). As emphasized by the authors of ref. |] the CP violating term 



in C reacts back on the skyrmion. We adopt their generalization of ansatz [2.6| and 
consider 

U c (r) = ex P (|[y|(/o(r)-5)J + A 8 / 8 (r)])C/ (r) (2.7) 

where 
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4 a + ml 



is the further vacuum shift required in order to prevent the creation of Goldstone 
bosons from the vacuum.^ The ansatz in eq. 277 remains spherically symmetric under 
combined isospin and spatial rotations. Neglecting isospin violation one can take 

,2 



Mi = M 2 



mz 



Mi 



2m 2 K — ml 



3a/N c 



2m\ 



(2i 



so that \g ~ Since it does not affect any issues of principle, we neglect 

77-77' mixing and corrections of order O(Mifa) in order to simplify the discussion. 
Taking these effects into account would only produce small numerical corrections 
of our results. The CP violating (CPV) part of the lagrangian, which contains 
contributions from both the mass and anomaly terms,§ provides a source for /o( r ) 
and f s (r) : 



5C 



CPV 



l-cosF)(/ 8 + v / 2/ ) + O( 



Mi, 



(2.9) 



2^3 v ' a 

The function F(r) is unchanged to order 0(d) by the presence of this perturbation. 

Neglecting terms of order 0(rn\jrn 2 ^) one can use this source to solve for f and 
fs in terms of F(r). The result isl 



fs(r) 
Mr) 



Xefn 

275 

2V3 



J3 / -m|r-r' 

— ^ -[l-cosF(r')] 

4tt |r-r'| 1 K n 

d 3 ?-' g— m'\r— r'| 

47r |r — r'l 



(2.10) 



[1 -cosF(r')] 



where we have neglected 77-77' mixing, and have written m = m n and m' = m rj i . The 
angular integration in eqs. |2.10| gives 



Mr) 
Mr) 



2\flm 
2^3- 



where 



/ dr'G m (r,r')[l -cosF(r')] 
Jo 

dr'G m ,(r, r') [1 - cos F(r')] 
e sinn(mr) r < r 



:2.iii 



m' Jo 







< 













G m (r,r') 

sinh (mr') r > r 

The contributions of /s( r ) an d /o( r ) will be important in what follows. 



(2.12) 



3. Meson-Nucleon Scattering in the Skyrme Model 

Here we generalize the analysis of ref. [3] to the CP- violating, three-flavour model 
described in the previous section. Denote the meson field (its quantum nature 
signaled by the absence of the subscript c) by 

U (r, t) = X a U a a = Q, 1,...,8. (3.1) 
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It is convenient to write 

U a = s a + ip a (3.2) 

where s a , and p a are real fields. The non-linear a model constraint that U is a U(3) 
field, 

U ] U = 1 , (3.3) 

gives nine constraint equations. Therefore we can regard the s a fields as dependent 
variables expressible in terms of the dynamical fields p a . (These nine independent 
real fields embody the pseudoscalar nonet.) 

Substituting eq. |3j] into the lagrangian (eq. $A\ j we obtain (exhibiting all time 
derivatives explicitly) 

C = -p a K aP p l3 + W a p a -V(p,Vp) (3.4) 

where a, (3 = 0,1,. ..,8 and W a is obtained from the Wess-Zumino term. (This is 
the analogue of eq. 3.3 of ref. ||) In the weak field limit 

K af3 = -^&a(3 + ■■■ (3.5) 

Beyond this, the detailed forms of K a/3 , W a , and V(p, Vp) will not be needed. The 
canonical momenta are 

6 L 

Tl a = — = K a pp p + W a (3.6) 

OPa 

so that 

p a = {K- l ) af3 {Up-Wp) . (3.7) 

Consequently, the Hamiltonian is, up to terms of 0(h 2 ) since we have ignored the 
problem of operator ordering, 

n= l -J d 3 x(n - w)^- 1 )^ - w)/3 + v(p, Vp) (3.8) 

and the canonical commutation relations are 

[p a (x) 1 U (3 (y)]5(x -y°)=t5 a ^ 4 \x-y) . (3.9) 

Due to the derivative couplings K aj 3 ^ -f&a.p and 7^ int 7^ —Ant; as is evident from 
eq. |3.4| - |3.9| . This is the essential reason that the classical argument for the absence 
of Yukawa couplings described in the introduction fails. 

The skyrmion soliton is the static but rotated solution (i.e. , U s = AUcA' 1 ) to 
the equations of motion obtained from eq. |3.4|. (The superscript s will denote the 



static, classical part of a quantity in the following.) In terms of the field of eq. \2.1\ . 
this solution is given by 

p s a = hmTr(X a AU c A- 1 ) (3.10) 
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where A e SU(3) rotates the ansatz U c to an arbitrary direction in SU(3) space and 
for the static solution we neglect any time dependence of A. 

In order to obtain the S-matrix for the meson-nucleon elastic scattering, we 
will analyze the LSZ reduction formula. To carry this out we must identify the 
interpolating field $ a . In fact, it is proportional to the quantum field p a in eq. |3.4| 
(not the classical p s a ). Since lim U c (r) = I, we can pick out the proportionality 
constant from the non-linear meson fields which appear in eqs. [L2] & |1.3| by using 
the fact that lim U(r,t) = U 2 d>: 

r—>oo Y 

+ + (3.11) 

so that 

$ a = yP« ■ (3.12) 

Therefore 

$ a =^ImTr[\ a U(r,t)] (3.13) 

where U(r,t) is the complete field U $AXJ C A~ X XJ $ that enters £(U). 

The reduction formula for our problem is given in terms of the source term J~(x), 
$ and $, as follows 

S fi = 5 fi + i 2 J dVdt' J dhdtft, {x')f h {x) 
x (B>(p')\T(J p (x')J a (x)) + 6{f - t)[<MA Ja(?)] - iu k> 6(t' - t)i<5>fi(x'), J a (x)]\B(p) 



(3-14) 

with fk(x) = e~ lk ' x I J (27r) 3 2c<jfc, where k ■ x = u^t — k ■ r. The source term of the 
meson fields is 

J a (x) = -[H, [H, $ a (x)]] + (-V 2 + ml)$ a {x) (3.15) 
with the commutator computed using eq. |3l| [3]^, and |3.12| . The result is 



+|(n - W)y(K~pK~£p - 2K^K- S ] p )(n - W)s (3- 16 ) 

+(n - w^K^KjiyVfu - Wsj) 

where for any quantity T and any index a, the notation T Q denotes ST/ 6p a . We have 
not kept track of operator orderings, so that eq. |3.16| is accurate only up to terms of 



0(h 2 ). We again see that the presence of derivative couplings in £ and H produces 



a different result in eq. 3.16 than that expected from classical considerations. 
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Similarly, using eq. |3.15| & |3.16| , we find that 

sj a (y) 



Mx),J a (y)]5(x°-y^ 



-%{K- X 



5(x° - y°) 



(3.17) 



(K- l ) Py \W y (x),J a (y)]S(x u -y l) ) 



with 



We also need 



[^(x),J a (y)]5(x°-y° 



Mx),J Q (y)}S(x°-y°) 



[$ (x) : J a (y)}5(x o -y°)= l 



2 SUJx) 



6(x" - y") 



(3.18) 



(3.19) 



The above equations are again accurate only to order h. 

In order to evaluate the Born terms from eq. 3.14| we need certain one-baryon 
matrix elements, which will be analyzed in terms of the 1/N C expansion. 

The crucial point is that in the matrix elements of operators taken between 
single-baryon states, one may replace the field operators with the classical fields 
evaluated with the skyrmion solution ( |3.10[) , i.e. , p a (x) — > P s a (x). This idea goes all 



the way back to reference |T2[ The reason is that the fluctuating parts vanish since 
the "in" and "out" meson states form a Fock space on the one-baryon subspace. 
As emphasized in ref. |] it is not necessary to explicitly seperate $ Q into static and 
fluctuating parts, since the one-baryon meson matrix elements will automatically 
pick out the static solution. 

In order to arrange the one-baryon matrix elements in powers of 1/N C , we need to 
order the various classical field quantities in terms of powers of iV c . By the classical 
equations of motion 



6p s 







Note that the ordering in 1/N C is the same for the CP conserving (CPC) and CP 
violating (CPV) terms in p s a , to leading order in CP violation, as can be seen from 



the following. From eq. |3.10|, the CP conserving term in p s a is 



(Pa)opc = -triXaArA- 1 ) • rsinF(r) 



so that 



I CPC 



0(1) 



(3.21) 
(3.22) 
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in the 1/N C expansion. Then, 

Pa 



(Pa)cPC + (Pa)cPV 

0(l)[l+X e /ml} 



(3.23) 



0(1)[1 + 



where we have used eqs. [O and Therefore, the l/iV c estimates that follow 

hold in both the CP violating and conserving cases. 



(3.24) 
(3.25) 



The estimates in eqs. |3.22| and |3.23| , and eq. |3J5| imply that 

(n.) ~ oil) 

Since 



one also has 



~ ^(iV- 1 ) (3.26) 

(jy% ~ o(i) 

for these classical quantities. With these ingredients we can now make 1/N C esti- 
mates of the one-baryon matrix elements. 



For the first term in eq. |3.14j we need 

(B\J a (x)\B>) = (B\(-V 2 + ml)® a + d?<S> a \B') 

= #{<5|(-V 2 +ml)(p s ) a \B>) + 0(1/N 2 C )} 



(3.27) 



which holds because all the terms in eq. |3.16j , evaluated with the classical solution, 
contribute at most to the order 0(1/ N^) correction. 
Similarly, from eq. |3.19| 



{B\[<$> p {x),J a {y)]5{x»-y°)\B') 



i f —{B\ 
2 x 1 



5 My) 
6U^{x) 



\B')5(x°-y°) 



0{f- A ) ~ 0{N e 



(3.28) 



and from eq. |3.18 
(B\[^(x),J a (y)}S(x -y°)\B') 



2 
0(1 



{-i(B\{K l ) 



Pi 



SJM 
Sp 7 (x) _ 



\B')5(x -y ) + O(l/N 2 c )} 
(3.29) 



8 



Therefore only eq. |3.27| and [3.29| contribute to the Born terms in the large N c limit. 
In addition, only equation 3.27 appears in the pole term, and so gives both the 
CP conserving and CP violating Yukawa couplings, while eq. |3.29| contributes to a 
non-pole contact term. 



4. The 7r-N Coupling Constants 



The classical, CP-violating solution is, from eqs. |2.7| , |3.12| , and |3.13 



(^)cpv = -^g {[V2(/o(r) - 5) + / 8 (r)](cosF(r) - 1) + 4/ 8 (r)} Tr(X a AX 8 A- 1 ) 

(4.1) 

Therefore, to leading order in 1/N C , 



g nNN (k) = J-^ e -^(B\(J a ) CPY \B>) 



(27T) 
d 3 



(4.2) 



(2tt) s 

If the physical (mass-shell) coupling 

VnNN = 9nNN{ k )W 2 =-ml (4.3) 

is to have a contribution from eq. |47|, then {$ s a )cpv must have a pion tail (z.e. it 
must behave as exp(— m n r)/r as r — ■+ oo). Since 

F(r) — > C- as r — ► oo (4.4) 

r 

it is clear that the cosF(r) factor does not have a pion tail. Further, eq. |2.10| shows 
that /s( r ) an d /o( r ) have (momentum space) poles at m 2 and mjL respectively, 
rather than at m%, so that the required pion tail is entirely absent from eq. |4.1| . By 
contrast, the CP conserving coupling is given by (the analysis is essentially the same 
as that already discussed in ref. |3|) 

2 k 2 + ?T7, 2 r 

9*nnQl) = 2 M *U | k | - J d 3 rjl (kr) sinF(r) (4.5) 

where ji(y) is the spherical Bessel function. Since sinF(r) (~ F(r) as r —>■ oo) does 
have the required pion tail, g n NN is no n- vanishing on-shell. 

Therefore we conclude that g nNN is suppressed relative to eq. [L6| by at least a 
factor of iV c -1 , while the behaviour of g n NN is correctly given by |1.6| .f| It is important 

t Although the order O(N^) contribution to H^nn does not vanish off-shell, it is suppressed by 



the factor m^/mz, as is evident from eq. 3.23. 
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to note that since g^NN is not suppressed by factors of N~ l , the stability of the linear 
fluctuations about the classical solution cannot be the reason for the suppression of 
'g-nNN- The correct explanation rests with the asymptotic behaviour of the skyrmion 
profile function eq. |2l| where one should compare the coefficient of As with that of 
r. 



5. Conclusions 

Our calculation shows that the Yukawa couplings in the Skyrme model involve 
issues more subtle than the vanishing of linear fluctuations about the classical so- 
lution. The CP conserving coupling g n n n is not suppressed, while the CP violating 
coupling is. Although this supports the final result of ref. |], we find that neglect of 
the pion contribution to D n in the large N c limit requires more analysis than indi- 
cated there before one can be assured that the direct contribution to D n dominates. 
The behaviour of the CP violating coupling constant g vNN is also interesting for 
other applications.il! 

The primary value of the Skyrme model is as a phenomenological approximation 
to the low-energy effective lagrangian of QCD; therefore, one hopes that it substan- 
tially embodies the physics of chiral symmetry and the large- N c properties of QCD. 
In fact, we have found that simple QCD estimates based on naive quark-line count- 



ing (i.e. estimates similar to those of ref. 14 but using the 6g FF interaction) agree 



with our Skyrme model results. In regard to this, it would have been peculiar if the 
Skyrme model, which does give the correct large N c behaviour for the multi-meson 
skyrmion, did not, in the end, give the correct Yukawa couplings as well. 
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